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Quantum Dots in a Strong Magnetic Field. 
Quasi-classical consideration 

A. MatulisQ 

Institute of Semiconductor Physics, Gostauto 11, 2600 Vilnius, Lithuania 
(Dated: February 1, 2008) 

The electron motion in rather strong magnetic fields (when only the lowest Landau level is popu- 
lated) is considered. In this case the electron kinetic energy is frozen out and the electrons are guided 
by slowly varied potential. Using the adiabatic procedure and expansion in magnetic length series 
the approximate description is developed. In zero order this approximation leads to the classical 
equations of motion describing the Larmor circle drift in the potential gradient. In the second order 
the special quantum mechanical description where the electron potential energy plays the role of 
the total Hamiltonian is constructed. Simple examples of a single and two electrons in the parabolic 
dot demonstrates that the proposed approximate description gives the main features of the electron 
system spectrum and the collective phenomena. 
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I. INTRODUCTION 

Quantum dots, or artificial atoms, have been a subject 
of intense theoretical and experimental research over the 
last few years jjj . The useful instrument in spectroscopy 
experiments is the magnetic field applied in perpendicu- 
lar to the quantum dot plane direction which enables to 
trace easily the dependence of the quantum dot proper- 
ties on various parameters. Moreover the strong mag- 
netic field reveals the quantization effects introducing 
into the electron system the favorable interplay between 
confining potential and Landau levels. 

Recently the main interest in quantum dots is re- 
lated to the electron-electron interaction and the collec- 
tive phenomena, such as the change of the ground state 
multiplicity, the electron density reconstruction, and the 
Wigner crystallization. The electron density reconstruc- 
tion in the finite electron systems was considered in [Q. 
Now it is known as Shamon-Wen edge — some of the elec- 
tron density ring around the finite system. Under certain 
circumstances the ring was reported to become unsta- 
ble H, and it breaks into separate lumps. Although the 
possibility to obtain the symmetry braking solutions was 
argued j| considering them as an artifacts of the approx- 
imate methods used, the exact calculations of the elec- 
tron correlation function Q undoubtedly indicates that 
the Wigner crystallization occurs at rather large electron- 
electron interaction. The presented in || electron den- 
sity plots show that the strong magnetic field facilitates 
the electron density edge reconstruction leading to the 
Wigner crystallization. 

Meanwhile the minimization of the system potential 
presented in JtJ shows that the Wigner crystallization in 
quantum dots can be successfully considered by classical 
mechanics. The fact that the strong magnetic field facil- 
itates the Wigner crystallization enables to suppose that 
the electron system behavior in very strong magnetic can 



be described by classical or quasi-classical methods. The 
purpose of the present paper is to show how such methods 
could be developed. The paper is organized as follows. 
After the formulation of the model in the next Section, 
in Sections 3 and 4 the main instrument — fast and slow 
variables are introduced. Then in Section 5 the adiabatic 
procedure is discussed and the slow motion Schrodinger 
equation is considered. In Section 6 the classical equa- 
tions for the limit case of strong magnetic field are de- 
rived, and in the next two Sections the illustrations of 
the simplified quantum mechanical description are given. 
In Appendix A the details of the adiabatic procedure are 
presented, and in Appendix B the transformation back 
to the initial coordinates is discussed. 



II. MODEL 



We consider the Schrodinger equation 



8 



ih—^f = H T ^ 

dt 



with the Hamiltonian 



Ht = ^-{p+-A(r)Y 
2m I c > 



V(r) 



(1) 



(2) 



describing the motion of 2D electrons in the strong per- 
pendicular homogeneous magnetic field and slowly vary- 
ing potential V(r). For the sake of simplicity the main 
equations will be derived for a single electron as the gen- 
eralization for the system of many electrons is trivial. It 
will be presented at the end of derivation. 

Choosing the symmetric gauge A = [B x r] we write 
down the main part of the Hamiltonian as follows: 



Ho 




eB 

Py + —X 



(3) 
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We shall consider it as a largest one treating the remain- 
ing potential V(R) as a small perturbation. 
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III. LANDAU LEVELS 

As in the standard perturbation technique we have to 
start with the zero order problem and solve the following 
stationary Schrodinger equation: 



{H - e}^j = 0. 



(4) 



The solution of it is known as Landau levels. The most 
simple way to obtain it is to introduce the new variables 



Pa 



1 

We 



-y, v 



-Py 



1 

2f F 



(5) 



where Ib = \JchjeB is the magnetic length. Using the 
new variables Hamiltonian (0) can be rewritten as 



TlLOr 



(6) 



where lo c = eB/mc is the cyclotron frequency, and the 
new variables obey the following commutation rule: 



(7) 



The zero order Hamiltonian reminds the Hamiltonian 
of the harmonic oscillator, and it is evident that it has 
the equidistant discrete spectrum which as it was already 
mentioned is called the Landau levels. 

We shall consider the case of very strong magnetic field 
when the electrons are in the lowest Landau level. Our 
task is to reveal how the slowly varying additional po- 
tential V(r) (as compared with the magnetic length Is) 
changes their behavior. 



IV. SLOW VARIABLES 

We shall treat the variables introduced in the previous 
Section as fast variables because they are included into 
the main part of the Hamiltonian. But as we are going 
to solve the 2D problem they are not sufficient to treat 
initial Schrodinger equation (Q). We have to introduce 
two more variables. We shall do that in the following 
way: 



1 1% 1 1% 



(8) 



We chose them in such way in order to have the most 
simple commutation relations, namely, 



[Z,X] = [Z,Y] = [r),X] = [r,,Y]=0, 



and 



[Y,X] 



-il 2 B . 



We shall consider those variables as slow ones. 
Now substituting the initial variables 

x = X + l B r), y = Y -l B £ 



(9) 



(10) 



(11) 



into Hamiltonian ■ 
sion: 



we arrive at the following expres- 



7 1 2 ) + V(X + I b V,Y~IbO- (12) 



So, we divided the Hamiltonian into two parts. The first 
largest one describing the motion of the electron in the 
homogeneous magnetic field depends on the fast variables 
only, while the other one — the slowly varying potential 
- depends on both fast and slow variables. Thus, we 
see that the slow and fast variables can not be separated 
exactly, but the presence of the small parameter (namely, 
the ratio of the magnetic length I B and the characteristic 
potential variation length lo ~ |U/VU|) enables us to 
separate them approximately by means of some adiabatic 
procedure. 



V. ADIABATIC PROCEDURE 

Now we are going to develop some adiabatic procedure 
and apply it for considering the Schrodinger equation (Q). 
For this purpose we shall make the following steps: 

• we expand the potential into ^-powers: 

V = V(X,Y) + IbvVx(X,Y) 

- I b £Vy(X,Y) + ---; (13) 

• divide the Hamiltonian into two parts: 

H = Hf + H s , (14) 

+ I b t]Vx(X,Y)+IbZVy(X,Y) + ---, (15) 
H s = V(X,Y); (16) 

• present the wave function as the product of its fast 
and slow parts: 



9 = il,(ri\X,Y)$(X), 



(17) 



• and use the following equation for the fast wave 
function part: 



{H f -E(X,Y)}iP(r,\X,Y)=0. 



(18) 



Actually it is the standard adiabatic procedure which 
has to lead to the Schrodinger equation for the slow elec- 
tron motion 

d 

ih— ®(X) = H<f>{X), (19) 

with the effective slow motion Hamiltonian 

H = V(X,Y) +E(X,Y). (20) 

However, there are some peculiarities caused by the 
fact that according to Eqs. (|?],[l0|) neither fast nor slow 
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variables commute each with other. That is why both 
wave function parts in Eq. ([l?]) depend only on a single 
variable (either r\ or X), while the other one has to be 
treated as an operator (£ = —id/drj, Y = —Igd/dX). 
Consequently, X and Y variables entering the fast wave 
function part %j){rf\X, Y) and the corresponding eigen- 
value E(X, Y) can not be treated as parameters (what is 
done in the standard adiabatic procedure), but should be 
considered as the operators acting on the slow wave func- 
tion part. This makes the adiabatic procedure a little bit 
tricky and cumbersome. Nevertheless due to the presence 
of the small parameter I b /Iq it can be performed. The de- 
tails of this derivation are presented in the Appendix [A|. 
Restricting the consideration up to the l B order we shall 
use the following slow motion Hamiltonian: 



the operators. Thus, we write 



dt 



;2 

l B 



- n [H,X] = ^ [ H,X] 



{H,X} = 



B 

c dV 
'eBdY' 



dt' 



-Y 



c dV 
^BdX' 



(24) 
(25) 



Note in the Heisenberg equations of motion the Plank 
constant fi is used (in spite of the fact that commutator 
of the variables is proportional to the magnetic length 
squared), because it has to be in agreement with the slow 
motion Schrodinger equation (19). Those two equations 
of motion can be rewritten as a single vector equation 



H = V^(R)+ l -fV 2 V^(R) 



(21) 



The superscript ^ indicates that the expression should 
be symmetrized in respect of the permutation of the slow 
variables X and Y which as we know already do not 
commute each with other. 

The above adiabatic procedure can be easily general- 
ized for the case of many electron system. As the slow 
motion different electron coordinates Ri commute each 
with other this generalization reduces to inserting the 
proper summations into obtained slow motion Hamilto- 
nian and replacing it by the following expression: 



N 



H = vWiRuRz, • • + f E V?^(fli, Jfe, • • •)• 

i=i 

(22) 

Now we are going to consider some simple examples in 
order to illustrate the application of the proposed simpli- 
fied description of the motion of electrons in the case of 
strong magnetic fields. Let us start with the zero order 
(Ib = 0) approximation. 



R= -\e z 



V]V(R) 



(26) 



where the symbol e z stands for the unit vector perpen- 
dicular to the electron motion plane z = 0. It is well 
known equation in plasma physics, and it describes the 
Larmor circle (the rotating electron in a strong magnetic 
field) drift caused by the gradient of applied additional 
potential. 

Thus, we see that system of 2D electrons in the very 
strong magnetic field (in the conditions of the fractal 
Hall effect, when only the part of the lowest Landau level 
is populated) demonstrates the classical behavior. This 
classical behavior is rather tricky. They do not behave 
as electrons. They behave as a system of classical gyro- 
scopes. 

Now let us go back and take the l 2 B order terms into 
account. In this case the quantum mechanical correction 
should take place, and we have to obtain something like 
quasi-classical description. In order to understand the 
main features of such quasi-classical motion let us take 
the most simple example of the parabolic dot with one 
and two electrons. 



VI. CLASSICAL EQUATIONS OF MOTION 



VII. SINGLE ELECTRON IN A PARABOLIC 
DOT 



In zero order approximation we shall take into account 
only the first term in Hamiltonian ( |2~]] ) and neglect the 
commutator ( |l0| ) between X and Y coordinates. We 
know that neglecting the commutators we have to arrive 
to the classical mechanics. But one has to remember that 
it is not correct just to neglect the commutators. It is 
necessary to replace them by the corresponding Poisson 
brackets according to the following rule (note we inserted 
l B instead of %): 



/2 
1 B 



[A,B] 



{A,B} 



dA dB 
dy dx 



dA dB 
dx dy 



(23) 



The most simple way to obtain the classical equations of 
motion is to use the Heisenberg equations of motion for 



In order to check the correctness of the above described 
method let us start with trivial problem of a single elec- 
tron in a parabolic dot. In this case we have the following 
potential: 



V(r) 



J 2 



(27) 



with the frequency ujq characterizing the strength of the 
confining potential, and according to Eq. ( pl| ) the follow- 
ing slow motion Hamiltonian: 



H 



J JL {X 2 



Y 2 )- 
d 2 

I 4 

lB dx 2 



1 



2/2 

muJ l B 



■X* 



U 2 



(28) 
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FIG. 1: Electron spectrum in a parabolic dot: solid curves 
— the exact result according dashed curves — the slow 
motion approximation (]2Sl). 



This the well known Hamiltonian of the harmonic oscilla- 
tor. Its eigenvalues and the corresponding eigenfunctions 
are 



(n + l), (29) 
- = e- x2 / 2l *H n (X/l B ) (30) 



where the parameter 7 = lo c /^q characterizes the relative 
strength of the magnetic field, and the symbol H n stands 
for the Hermit polynomial. 

In order to evaluate the approximation obtained by 
solving the slow motion Schrodinger equation let us com- 
pare it with the exact Fock-Darvin result which is 

|(2n + |m| + l)V!+7 2 /4+ (to - 1)7/2} 

(31) 

where orbital quantum number to is integer, and radial 
quantum number n is integer and nonnegative. This ex- 
act result together with the approximate one (^) are 
shown in Fig. [I] by solid and dashed curves, correspond- 
ingly. We see that in the asymptotic region 7 — > 00 
(shown by the dotted rectangular) the approximate re- 
sult is rather close to the rotational levels belonging to 
the lowest Landau level. Moreover, we may expect the 
quantitative agreement already at 7 > 2 values. It is 
interesting to inspect how the wave function and the cor- 
responding electron density looks like. However, we have 
to remember that eigenfunction ( p0| ) is not the electron 
wave function itself but it is its slow motion part only. 
In order to obtain the electron wave function according 
to Eq. ( |l7| ) we have to multiply it by fast motion part. 
Next we have to go back to the initial variables (pd|). It 



can be done using some integral transformation which is 
described in Appendix Using transformation kernel 
( p36| ) and restricting our consideration by the lowest fast 
wave function approximation flAl6|) we write down the 
total electron wave function in initial x, y variables 



*n(x,v)= dr) dX(x 1 y\ Vl X)Mv)^n(X) 

J — OO J — OO 

= t= \ dn dXHJX/l B ) 

2ttI b V2^.J-co J -00 ' ' 

. e -r,*/2 + iy(X-l B r,)/2ll-X*/2ll + ^ _ x) 

(ixy-x 2 )/2l 2 b roc 

drje-'i +{x - ty ^ ,lB H n {x/l B - 7]). 

(32) 



2irl B V2 n nl 



Fortunately the integral can be calculated by analytical 
means. Using the standard integrals with Hermit polyno- 
mials H we obtain the following expression for the total 
electron wave function 



*«(a?,v) = 



1 



l B \ / 2 n + 1 n\n 



i mv {r/l B ) n e r2 / 4l B, (33) 



and the corresponding electron density in the n eigenstate 
Pn{r)~{r/l B ) 2n c- r2 / 2l2 B. (34) 

We see that in the case of large n values (in the quasi- 
classical case) the electrons are mainly located on the 
ring. Equating to zero the derivative of the above density 
expression we obtain the radius of this ring. It is 



r = l B V2n. (35) 
Now inserting the n value expressed from Eq. (|29|) we get 



ro = h 



2-fE 

hujQ 



2E 



(36) 



what exactly corresponds to the classical potential energy 
E = V(ro) = ttlljqTq j2 of the rotating electron drifting in 
the confining potential along the circle with the radius ro- 
The single difference of quasi-classical electron behavior 
from the classical one is that now according to Eq. (|34|) 
it moves not along the thin trajectory, but it is spread 
over the ring with the thickness of order / B . 



VIII. TWO ELECTRONS IN A DOT 

The other good for us example is the two electrons in 
a parabolic dot because there is the exact solution which 
can be compared with our approximate results (see, for 
instance, ||). In this case the behavior of electrons is 
described by the following potential: 



V(n,r 2 ) = ^{r\+rl} + 



n - r 2 \ 



(37) 
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Let us introduce the center of mass and relative motion 
coordinates. We shall do it in a non standard way in 
order not to spoil the commutations rules for the fast 
and slow variables which we already used. Namely, we 
use the following definition: 

r c = 4=(n + r 2 ), r r = i(r 1 -r 2 ). (38) 

It leads to the separation of variables as the potential can 
be presented as a sum of two terms 



V(r 1 ,r 2 ) = V c (r c ) + V I (r I ). 
The potential for the center of mass motion 



(39) 



Consequently, the relative slow motion eigenvalue reads 



e { ; ] = huso 



Av/7 



7 y/2{2n+ 1) 



4(2n + l) 



where dimensionless parameter of electron-electron inter- 
action A = la/as is the ratio of the characteristic confin- 
ing potential length Iq = y^h/mujQ and the Bohr radius 
as = h 2 /me 2 . 

Now adding together eigenvalue (|2^) for center of mass 
motion, the relative motion eigenvalue (f45|), and one 
more term huj c for lowest Landau level energy we ob- 
tain the final result for the two electron eigenvalue in the 
parabolic dot in the slow motion approximation 



muj. 



0^2 



(40) 



V c (r c ) 



exactly coincides with the single electron potential (|27^ 
which was already considered in the previous Section. 
Consequently, the eigenvalue and eigenfunction of the 
center of mass motion coincide with those given by 
Eqs. ( |29| , p0| ). Note that now the capitals X and Y have 
to be replaced by the slow center of mass motion coor- 
dinates X c and Y c . Performing the same procedure as 
in the previous Section we shall arrive at center of mass 
motion density given by Eq. ( |34l) with the coordinate r 
replaced by the center of mass coordinate r c . 
The relative motion potential is given as follows: 



V t (r r ) = 



U 2 



(41) 



According to Eq. (^l|) it leads to the following slow rela- 
tive motion Hamiltonian: 



771Cl>q / 2 l B 



(42) 



The symbol R 2 of course has to be replaced by the oper- 
ator 



R r — ► — I 



d 2 



B dX} 



X 2 . 



(43) 



Consistently the slow motion Schrodinger equation with 
Hamiltonian (^) can be solved by means of Fourier 
transformation technique presented in Appendix [A|. But 
in this simple case of two electrons one can find the eigen- 
values of the above slow motion Hamiltonian rather eas- 
ily paying attention to the fact that the eigenfunctions 
of R 2 operator diagonalize Hamiltonian (^) as well. We 
know the eigenvalues and eigenfunctions of operator R 2 
already. They are given by Eqs. (^9|,[50|). Thus, in order 
to obtain the eigenvalues of Hamiltonian (|4|) we have to 
make the following replacement just in the above Hamil- 
tonian: 



R. 



Z|(2n + 1). 



(44) 



En.i 



N + n + 2 



7 



Av/7 



V2(2n+ 1) 



1 



4(2n + l) 



(46) 



The dimensionless eigenvalue (in Hujq units) dependencies 
on the relative magnetic field strengths (on parameter 
7 = lu c /luq) for the case of N — and several n values are 
shown in Fig. ||a. In Fig. |^b these eigenvalues are com- 
pared with the exact solution taken from Merkt paper [|| . 
We see that when n — > oo and 7 — > 00 (namely, in the 
asymptotic region) the approximate consideration is in 
good agreement with the exact one. Moreover, the quasi- 
classical treatment describes correctly the main features 
of the electron behavior in strong magnetic fields, namely, 
the increment of the angular momentum (the quantum 
number n plays its role) with the increment of the mag- 
netic field strength. Indeed, minimizing relative motion 
eigenvalue ( pl5| ) in respect the magnetic field 7 in the case 
of large quantum numbers n we obtain the ground state 
orbital momentum 



no = (A/4) 2 / 3 7 



(47) 



proportional to the magnetic field what agrees with the 
quantum mechanical result obtained in flCj] . 

Let us also look at the electron density given by the 
quasi-classical approximation. First, we notice that ac- 
cording what was said above the eigenfunction of Hamil- 
tonian ([l2]) coincides with the eigenfunction of the oper- 
ator R 2 . Thus, it coincides with wave function ( |30| ) with 
the coordinate X replaced by the relative motion coor- 
dinate X r . So, perfor min g the same transformation as it 
was done in Section VII we obtain the relative motion 



density given by expression (|34|). 

Now taking into account that N = corresponds to 
the ground two electon state we can write down the two 
electron distribution function in the following form: 



Pn(ri,r 2 ) 



2n p -r r 2 /2/ 2 f 



r 2 



{tt/Ib) 



|2» p -(r 2 1 +rl)/2l 2 l 



(48) 
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FIG. 2: Spectrum of two electron system in a parabolic dot: 
a — slow motion approximation results; b — the comparison 
of approximate results (dashed curves) with the exact results 
(solid curves) taken from Q. 



In Fig. ^ the above function is plotted as a function of 
the first electron coordinate with the other electron co- 
ordinate fixed at the point corresponding to its classical 
equilibrium position. This point is indicated by solid 
dot. Actually the plot represents so called pair correla- 
tion function. The plot in Fig. [|a corresponds to n = 5 
and the plot in Fig. ||b — to n = 2. We see that for 
larger n the pair correlation function demonstrates the 
peak in the opposite to the fixed electron position what 
corresponds to the Wigner crystallization of this simple 
two electron system in the strong magnetic field. 

When the magnetic field strength decreases (what cor- 
responds to the ground state with the smaller angular 
momentum value, say, n = 2 as it is shown in Fig. [|b) 
the Wigner crystal starts to melt — the pair correla- 
tion function transforms itself from the peak into ring. 
Note, it demonstrates the fact that the angular melting 



FIG. 3: Pair-correlation function for various orbital momenta. 



precedes the radial one what agrees with quantum me- 
chanical result in |0| obtained for the case without the 
magnetic field. 



IX. CONCLUSIONS 

Let us summarize shortly what was said about the be- 
havior of electrons in the case of the strong magnetic 
field. 

In the asymptotic region of the strong magnetic fields, 
when all electrons occupy the lowest Landau level only, 
their kinetic energy is frozen out, and their behavior is 
guided by weakly varying (characterized by some charac- 
teristic length Iq) additional potential. Applying special 
fast and slow motion variables, the adiabatic procedure 
and the expansion in Ib/Iq powers some simplified ap- 
proximate description can be developed. 

In zero order approximation one get the classical equa- 
tions which describe the electrons as a system of gyro- 
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scopes. Those equations actually are the equations for 
the Larmor circle drift in the gradient of applied poten- 
tial. 

Applying the expansion up to (Ib/Iq) 2 order one ob- 
tain the self consistent equation set which coincides with 
the Schrodinger equation where the role of canonical vari- 
ables play two cartesian slow motion 2D electron coordi- 
nates X and Y with the commutator proportional to the 
magnetic length squared (instead of being proportional 
to the Plank constant as it is in the standard Schrodinger 
equation). 

Two simple examples of a single and two electron in 
a parabolic dot demonstrate the accuracy and main fea- 
tures of proposed approximate description. 

The approximate eigenvalues coincide with the exact 
ones in the asymptotic region when 7 = lo c /ujq — > 00, and 
even in the intermediate region (7 > 2) one can expect 
rather good semi-quantitative description. 

The electron wave functions can be obtained as prod- 
uct of the fast wave function part corresponding to the 
lowest Landau level and slow wave function part obtained 
by the above specific slow motion Schrodinger equation. 
After the transformation to the original variables the 
wave function obtained in this way describes correctly 
the quantum mechanical Ib/Iq order correction to the 
classical electron motion. 

The two electron in a dot example shows that the pro- 
posed approximate consideration describes such collec- 
tive phenomena as the Wigner crystallization, the change 
of the angular momentum in the ground state when the 
magnetic field strength increases, the phenomena of an- 
gular and radial melting of the Wigner crystal. We hope 
that this approximate method can be useful for the con- 
sideration of more sophisticated many electron systems 
when the straightforward solution of quantum mechani- 
cal equation meets computational difficulties. 
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APPENDIX A: SLOW MOTION HAMILTONIAN 

As it was mentioned in Section [V| performing the adi- 
abatic procedure steps we have to pay attention to the 
fact that the variables X and Y do not commute each 
with other. That is why instead of using the straight- 
forward expansion (|l|) we apply the following Fourier 
transformation: 



V(x,y) = 




. j e ifc* e i?j, + e ^e lfcr }u(fc, q), (Al) 

dx / dye- ikx e- iqv V(x,y). (A2) 
-00 J — 00 



These two expressions can be considered as a definition 
of the operator function V(x,y). Thus in the first ex- 
pression the symbols x and y will be considered as the 
operators, while in the second one x and y are just the 
dummy integration variables. The main advantage of 
such potential representation is that the operators x and 
y are moved from the general potential function V(r) to 
more simple exponent functions. Although the old x and 
y variables commute we used the symmetric exponent 
product which will be necessary in further derivation. 

Now we substitute variables (|ll]) into exponents and 
expand them into £, 77-powers: 



e ik(X+l B ri) e iq(Y-l B £,) _ Q ikX Q iqY Q il B kj] B q(, 

= e ^e^{l + il B kr, - \}%k 2 rf}{l - U B qt 

= e lkX e^ Y {l + il B (kr, - qO - \l%k 2 rf - + l 2 B kq^}, (A3) 

e iq(Y~l B £) e ik(X+l B ri) _ & iqY Q ikX & -il B q^^il B kr) 

= e iqY e ikX {l + il B (k V - q£) - \l%kW - ^Wi* + l 2 B kq^}. (A4) 

Taking the evident equality 

a V li + bl;T) = ^a{rt + + + \b^V + rt-i) = \{a + b)(£r) + t£) + % -{a - b) (A5) 
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into account we write down the following expansion of the symmetric product of exponents: 

e ik(x+i B r,) e i<i(y-lB$) + e iq(Y-lBt) e ik(x-iBv) = f e ikX e iqY + e lqY e ikx } 

x {l + il B (kr, - qO - h%(k v - qO 2 } + ^|kz{e**e* y - e^e j «} 



JkXiqY 

e e v 



s ). 



e ikX e iqY 



(A) 



(A6) 



Note how the symmetric and antisymmetric exponent products and function 77, fc, g) are defined. 

Inserting the above expansion into Fourier transformation (Al), changing the parameters k and q by the operators 
id/dx and id/dy acting on exponents, and performing the integration by parts we arrive at the following potential 
expansion: 

fOG />00 



V(X + l BV ,Y-l B Z)= dx dy — ^-V{x,y) 

J—oo J—oo J — oo ^ J — oo ^ 

x \ L{i,ri,k,q) 



e tk(X-x) e tq(Y-y) 



(S) i 

+ 2 



l 2 B kq 



e ik(X-x) e 2q(Y-y) 



(A) 



OO POO 



,Lr j ,1/1 j ^- j J-V( x ,y)x I L(£,r), id/dx, id/dy) 



oo J — oo 



e ik{X-x) e iq(Y-y) 



(S) 



ill d 2 

2 dxdy 



e ik(X-x) e iq(Y-y) 



(A) 



dxldyl^f^ 



2tt / 2?r 



e ik(X-x) e iq(Y-y) 



s ). 



L(£,r),-id/dx, -id/dy) 



it 2 , 



e ik(X-x) e iq(Y-y) 



(A) d 2 



dxdy 



■V(x,y) 



dx I dy I — [ — — < 



2tt / 2?r 



e ik{X-x) e iq{Y-y) 



(S) 



I \ dx dy/ 2 V dx dy) J 



'dy 



e ik(X-x) e iq{Y-y) 



(A) d 2 



dxdy 



■V(x,y). 



(A7) 



Actually it is the definition of the operator function expansion which we have to use instead of expression (|13|). It 
can be rewritten in more simple way if we use the following operator function definition: 



F^ XY \x,Y)=Jdx JdyJ^J ^ x -^ Y -^F(x,y). 



(A8) 



It defines the function with ordered operators — all op- 
erators X stand on the left side of the operators Y in all 
terms of its Taylor expansion, or Fourier transform. 

Defining the symmetric and antisymmetric operator 
function as 

F^(X,Y) = 1{F& Y \X,Y) +F( YX HX,Y)}, (A9) 
F^ (A, Y) = \{f (XY) (A, Y) - F( YX 1 (X, Y) } (A10) 

we rewrite the potential expansion in the following formal 



simple form: 




(All) 



Now we are ready to perform the next step of our adia- 
batic procedure, namely, to insert the obtained potential 
expansion into fast Hamiltonian ( p"5| ) and solve the fast 
eigenvalue problem (|l^) . It can be easily performed using 
the standard perturbation technique. 
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Indeed, using the modified fast Hamiltonian 



Hf 


— Ho + Hi + H 2l 


(A12) 


Ho 


= ^« 2 + ^), 


(A13) 


Hi 




(A14) 


Hi 


- 4(4-4)'^ 


(A15) 



we obtain the zero order eigenvalue and function 

E Q = ^, Mv)=K- 1/4 e- n2/2 - (A16) 



Next, due to the zero order function symmetry we get 
Ei = 0, and solving the first order equation 

{H a - E a } = -H^o = -1 B (VP ~ ivj S) ) V1 p (A17) 
we define the first order correction to the wave function 



Mv\x,Y) = - lB ^]- iV ^\ ^. 

nuj r 



(A18) 



Then from the second order equation one can easily get 
the following second order eigenvalue correction: 



/oo 
dvMv)H 2 Mv) 
-00 



4 I 



y(S) , y(S) 
V XX ' V YY 



d v Mv)HiMv\x, Y) 

I 2 



(S)2 
Y 



}.(A19) 



Now having the fast motion problem eigenvalue calcu- 
lated we can proceed with the adiabatic procedure. For 
this purpose we present the total wave function as 

*(ri,X,t) = c-^^Mri+MvlXiY)}^^), (A20) 

insert it into Eq. (Q) , and obtain the following expression: 

• {ip (r))+Mv\X,Y)}$(X,t)=Q. (A21) 

Now multiplying the above equation by function ipoirj) 
from the left side, integrating it over all — 00 < ?; < 00 
interval and taking into account the action of the fast 
Hamiltonian Hf on the fast wave function part we ar- 
rive at the slow motion equation (|l9|) with the effective 
Hamiltonian 

h = + l ^vi A j + l ^v 2 v {s) - -^L{w( s n 2 . 

2 XY 4 2huj c L ' 

(A22) 



The performed procedure is consistent at least with the 
accuracy up to l 2 B . That is why we shall omit the second 
and the last terms in ( A22J) as they are of order l B . In the 
last term this dependence appears due to the additional 
factor uj c in the denominator, while in the second term 
it is caused by slow variable commutator (|l^). Omitting 
these terms we arrive at the final slow motion Hamilto- 
nian (pl|). 



APPENDIX B: COORDINATE 
TRANSFORMATION 



According to the ideas of the quantum mechanics we 
can change the wave function variables by means of the 
following transformation: 



d v dX{x,y\ii,X)*(n,X) (Bl) 



where the transformation function (x,y\r], X) has to be 
chosen as the eigenfunction of operators x and y with the 
corresponding eigenvalues x and y. Namely, this trans- 
formation function has to obey the following equations: 



{x - x}(x,y\r),X) = 0, 

{y-y}(x,y\v, x ) = o, 



(B2) 
(B3) 



{X + l BV -x}{x,y\i 1 ,X) = 0, (B4) 
-il%—+il B —-yUx,y\r,,X) = 0. (B5) 



It can be checked straightforwardly that the following 
transformation function satisfies both equations: 



(x,y\rj,X) = -^=e lv{x - lB,l)/2l ^S(X + l B T 1 -x). (B6) 



The normalization factor is chosen in agreement with the 
condition 



00 POO 



dx / dy (x,y\r),X)(x,y\r) ,X ) 

3 J — OO 

= S(r ] -r ) >)5(X-X>). (B7) 
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